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Population III Gamma Ray Bursts 

P. Meszaros 1 and M.J. Rees 2 

ABSTRACT 

We discuss a model of Poynting-dominated gamma-ray bursts from the 
collapse of very massive first generation (pop. Ill) stars. From redshifts 
of order 20, the resulting relativistic jets would radiate in the hard X-ray 
range around 50 keV and above, followed after roughly a day by an external 
shock component peaking around a few keV. On the same timescales an 
inverse Compton component around 75 GeV may be expected, as well as a 
possible infra-red flash. The fluences of these components would be above 
the threshold for detectors such as Swift and Fermi, providing potentially 
valuable information on the formation and properties of what may be the 
first luminous objects and their black holes in the high redshift Universe. 

Subject headings: gamma-rays: bursts — cosmology — stars: population III 
— jets: magnetized - radiation mechanisms: non-thermal 



1. Introduction 

Population III stars are widely considered to consist mainly of 'very massive stars' 
(VMS) in the the range of hundreds of solar masses (Ohkubo et al.; Yoshida et al. 
2006). The VMS are expected to be very fast rotating, close to the break-up speed, 
and accretion leads to a mass upper limit which may be around 1O 3 M . Those in 
the 14OM £ M* £ 260M© range are expected to be subject to pair instability and 
explode as supernovae without leaving any compact remnant behind, while those above 
~ 26OM are expected to undergo a core collapse leading directly to a central black 
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hole (Heger and Woosley 2002), whose mass would itself be hundreds of stellar masses. 
Accretion onto such massive black holes could lead to a scaled-up collapsar gamma-ray 
burst (Heger et al. 2003; Komissarov and Barkov 2009). In this paper we discuss a 
specific scenario for pop. Ill VMS collapsars at redshifts of order z ~ 20, resulting in 
Poynting dominated relativistic jets which produce GRBs with characteristic radiation 
properties extending from soft X-rays to multi-GeV energies. 



2. A population III collapsar model 

We consider a pop. Ill star undergoing core collapse at a redshift z ~ 20, which 
leaves behind a black hole of mass Mh surrounded by an accretion disk or torus of 
mass Md- Knowledge about the progenitor structure and the collapse history is rather 
limited, an approximate but representative scenario having been outlined in Komissarov 
and Barkov (2009). This assumes a nominal VMS of mass M* = 10 3 M 3 M Q and radius 
i?* = 10 12 -R* 5 i2 cm rotating at half the break-up speed, which results in a disk of outer 
radius of Rd, disk mass M d and central black hole of mass M h . Given the uncertainties, 
for the purposes of estimates we can assume that the typical disk mass and the black 
hole mass are of the same order as the progenitor mass, M d ~ M h ~ M* ~ 1O 3 M 3 M . 

For such large BH masses the accretion torus density and temperature are too low 
for neutrino cooling to be important, and the accretion regime can be described through 
an advection dominate (ADAF) model, e.g. Narayan and Yi (1994), in which radiation 
pressure is dominant. The properties of the precursor VMS, and the flow dynamics 
after the collapse, are plainly uncertain; it is nonetheless helpful to parametrize the key 
numbers and scaling relations in terms of an undoubtedly oversimplified but specific 
model. For a VMS rotating at, say, half the break-up speed the disk outer radius will 
be Rd = i?*/4, and for a disk viscosity parameter a = 10 _1 a;_i the accretion time 
t d ~ (U/9a)(R 3 /GMy/ 2 is approximately 

t ac ~ 5 x ltfaZ\Rl l2 M- 112 s, (1) 

and the mean accretion rate is M ~ ~ 0.2a^iR~ 3 / 2 2 M^ 2 M Q s -1 . The disk inner ra- 
dius at the marginally bound orbit of a rotating Kerr black hole with rotation parameter 
a is 

n = {R s /2)h{a) ~ 3 x 10 8 M 3 cm, (2) 

where R s = (2GM h /c 2 ) and f\(a) = [2 - a + 2(1 - a) 1 / 2 } ~ 2.1 is estimated for a ~ 0.8. 
The inefficient neutrino cooling is insufficient to power a strong jet, but strong magnetic 
field build-up in the torus could lead to much stronger MHD jets. The disk mass density p 
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and gas pressure P in the ADAF regime provide an estimate of the disk poloidal magnetic 
field as B 2 = (8tiP/(3), where /3 = 10/?i is the magnetization parameter, leading to a 
Blandford-Znajek type (Blandford and Znajek 1977) Poynting jet luminosity 

where f 2 (a) = (a/2)(l + y/1 - a 2 ) and for 0.5 £ a £ 1 the product ff 2 f 2 ~ 1/4 
(Komissarov and Barkov 2009). For a jet solid angle Q = 10 _2 f2_2, as may be expected 
in a large star, and assuming that a fraction 77 = 10~ 1 ?7 1 of the luminosity falls in the 
X-ray band (see §5), from a redshift z ~ 20 and using current cosmological parameters, 
using L iso = L(47r/tt) a flux F = r]L iso / %ixr\ ~ lO-^^QzlR'^M^ 2 erg cm" 2 s" 1 
is expected, more than an order of magnitude above the Swift BAT sensitivity. In the 
next section we show that such objects could indeed have such X-ray luminosities, as 
well as other components at higher energies. 



3. Population III Poynting dominated GRBs 

Taking the magnetic luminosity of eq. (3) as representative for a pop. Ill collapsar 
at a redshift z ~ 20, we use this as the central engine underlying a scaled-up version 
of a Poynting-dominated GRB model discussed in Meszaros and Rees, 1997; henceforth 
MR97). That is, we assume that a purely MHD (Poynting dominated) jet emanates 
from the central black hole plus accretion torus system, which is initially devoid of 
baryons. At the base of the jet given by eq. (2) the initial jet bulk Lorentz factor is 
parametrized as ~ L/L w , where L is the magnetic luminosity (3) and L w = L e ± rj is 
the associated pair and photon luminosity (allowing for the possibility that > 1, as in 
pulsar wind models). Depending on uncertain details at the base of the jet, the initial T e 
could differ from unity, T e <^ 1. The jet magnetic luminosity L is related to the comoving 
transverse magnetic field B' e at the base of the jet through L ~ 4irrjc(B' 2 /8ir)Tj, where 
L = 5 x 10 52 L 52 . 7 erg s^ 1 = 5 x 10 52 R~% 2 M^ 2 erg s _1 can be taken as the isotropic 
equivalent luminosity for an observed within an angle 1/T of a jet axis whose final 
Lorentz factor T ^> 1. 

We review here the dynamics of baryon-free Poynting jets for the VMS collapsar 
case, since the numerical values differ from those of the normal stellar collapsar case. 
Near the base ri of the outflow the jet will generally become loaded with pairs which 
are in near thermal equilibrium with photons. At the base of the jet the transverse 
comoving magnetic field strength, the comoving pair temperature and the comoving 
pair density are B\ ~ 6 x lO 12 ^^^ 1 G, T[ ~ 3.7 x l^Ll^r^Tf 1 ' K 
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and n' e ~ 1.5 x lO^Ls^^ 1 * cm" 3 , where T £ ' ~ (L w / ' ATXcr 2 Y 2 a B f/ A and n£ ~ 
(0Aa B kT' 4 /kT£), with a# is the Stefan-Boltzmann constant, 77,8.5 = (l/3 x 10 8 M3 cm) 
is given by eq. (2) and n' ~ n' ± ~ with T' ± ~ T^. The optically thick pairs and photons 
will be frozen-in with the magnetic field, the whole behaving as a relativistic fluid. The 
magnetic geometry is that of an aligned rotator, the field spiraling out along the axis, 
and the pair plasma is frozen-in. Assuming for simplicity a constant opening angle jet 
with purely radial motion, the lab frame transverse magnetic field strength varies as 
B oc r -1 , and the comoving transverse magnetic field strength varies as B' oc B/Y. 
Along the jet the transverse field lines do not change polarity, which is less likely to lead 
to reconnection, especially for an essentially baryon-free outflow such as we assume here. 
This differs from models, e.g. Drenkhahn and Spruit (2002), which assume a baryon load 
and where reconnection plays a role in the dynamics of the expansion. In the absence of 
reconnection, the magnetically dominated comoving energy density e' oc r" 2 r~ 2 , while 
the comoving pair density n' oc T' 3 oc r~ 3 . The fluid has a bulk Lorentz factor Y in 
the lab frame and an internal (comoving) Lorentz factor 7' ~ s'/n', and from energy 
and entropy conservation, Y grows at the expense of 7', so Y oc 1/7' oc n'/e' oc r, for 
a jet of constant opening angle. Hence, initially Y ~ Y^r/rc), and the comoving pair 
equilibrium temperature varies as T' oc r" 1 . 

When the comoving temperature drops below m e c 2 the pairs start to recombine, but 
the linear acceleration continues beyond this point, until the comoving pair Thomson op- 
tical depth t' t has become less than unity, which occurs when the comoving temperature 
V oc r- 1 has dropped to kT' a ~ 0.04m e c 2 ~ 17 keV (e.g. Shemi and Piran, 1990). This 
occurs in our case at a radius r a where (r a /r^) = (T' t /T' a ) = (r a /r^) ~ 20L52 4 7 r^Tg ^YJ 3 ^ 4 , 
or 

r a * 6 x 10%^Y^ , 

r a ^ 2 x lo^L^r-^rf . (4) 

At this radius most of the pairs have already recombined, and the gas density consist 
mainly of photons, n' a ~ n' e (r a /re)~ 3 ~ 2 x 10 26 cm" 3 , while the remaining pair comoving 
density, from Saha's equation, is n'± a ~ 5 x 10 15 r^r^g 5 cm" 3 . The photon to pair ratio 
at this radius is ~ 4 x 10 10 , and pair annihilation practically ceases beyond this, so 
the maximum theoretical (inertial limit) Lorentz factor is Y in ~ (n^/n' ± )(kT^/m e c 2 ) ~ 
1.6 x 10 9 . Other effects, however, can set in before that, resulting in a lower terminal 
value. 

Above the annihilation radius r a , the comoving density of pairs n' = n'± providing 
the inertia has been drastically reduced, but the magnetic pressure or the comoving 
energy density e' oc B' 2 continues acting continuously. In the lab frame the transverse 
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component of the field varies as B oc r" 1 . Since T oc n'/e' oc n'/B' oc B 2 /n', the 
drop in n' implies that the gas must accelerate faster than the previous behavior of 
T oc r. The pair density is above the minimum below wich the plasma can't carry the 
MHD currents even if the stream velocity of the positive relative to negative charges is 
of order c (Goldreich and Juian 1969), which is essentially the same for an e + e~ or an 
e~p + plasma (note that the Blandford-Znajek mechanism exploits the analogy with the 
aligned pulsar case); thus the MHD regime remains valid. 

In the next three paragraphs we outline a possible scenario for the dynamics above 
the photospheric radius of a Poynting jet, which can lead to very high bulk Lorentz 
factors T. This high-r scenario is sensitive to model details, and must be considered 
speculative. For instance, such high Lorentz factors could lead to gradients which inval- 
idate the transverse field assumption; also, a small amount of entrained baryons could 
end up dominating the inertia. For simplicity we neglect such potential complications, 
noting that our final observational predictions are essentially idependendent of what 
happens in this high-r regime. Thus, we consider a highly idealized picture of what 
happens above the photospheric radius r a . While most of the photons escape freely, 
the pairs continue being scattered repeatedly by the much more numerous photons, and 
continue experiencing a drag for some distance r > r a . It is useful to compute this in a 
frame Tj ~ (r/r^)T^ where the photons are isotropic. In this frame the electrons, which 
are essentially cold in the comoving frame, are boosted to a Lorentz factor 7 ~ T/Ti, 
and the drag time is ta r ,i ~ m e c 2 / \u p h^<j T c~f) = (m e c 2 4:7ir 2 Tf/LaTT), where u p h is the 
radiation energy density. In the lab frame the drag time is Tj times longer, and the 
acceleration rate is obtained by setting the ratio of the lab frame Compton drag time 
tdr = (m e c 2 4iTr 2 / Lo"jT)(r /r^) 4 and the lab frame expansion time r/c equal to the ratio of 
the kinetic and the Poynting flux, n'm e c 2 T 2 /[{Bf/4:Ti){r^/r) 2 ]. The drag time is shorter 
than the annihilation time at r > r a so the remaining pairs are frozen in, and for a 
dominantly transverse magnetic field, the ratio of the drag to expansion time is oc r 5 /T 
while the ratio of the kinetic to Poynting flux is oc T, so for r > r a the flow accelerates 
as T oc r 5/2 . 

As the Lorentz factor continues increasing beyond r a the annihilation photons, 
whose isotropic frame energy is 0.12m e c 2 (r a /r), eventually are blueshifted in the jet 
frame to ^ m e c 2 , and their directions are randomized by scattering. Given that the 
compactness parameter is large, this results in copious pair production 77 — > e + e~ . This 
reduces the drag while increasing the inertia, leading to a mitigation of the acceleration 
rate beyond a radius r p where the Lorentz factor is r p , 

r p ~ 4.6 x lO^VJ^r,- 3 ' 4 cm, 

r p ^ 1.5 x lo^^r^rf- (5) 
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Beyond r p , pair formation will be self-limiting, resulting in a scattering optical depth 
t± ~ 1, the threshold condition being 0.12m e c 2 (r a /r)(r /Tj) ~ m e c 2 , and since Tj oc r 
this implies at r > r p an acceleration rate T oc r 2 . 

A basic uncertainty is that in addition to the jet's own annihilation photons there 
will be other photon sources which may also provide drag on the pairs. One such 
additional photon source is the accretion torus, but this radiation will be absorbed by 
the rest of the stellar envelope through which the jet is making its way, unless all of it is 
incorporated in the torus or is blown away. Also, the radius r p is near the outer envelope 
of a VMS of radius R* ~ 10 12 cm, which is a source of UV and soft X-ray photons. The 
drag on these external photons is likely to limit the Lorentz factor to values not much 
above T p . 

Irrespective of the final Lorentz reached in the previous high-r phase, outside the 
stellar radius R* the jet is expected to propagate through a stellar wind, whose details are 
poorly known for pop. Ill objects, and further out the jet will encounter the interstellar 
medium of the minihalo or protogalaxy hosting the VMS. The jet will shock and sweep 
up the external medium, pushing it ahead of itself. While the jet continues to be fed 
by accretion, over a time of order t ac given by eq. (1), after an initial brief transient 
the shocked jet head will continue to advance at a decelerating pace into the external 
medium. The Lorentz factor of the jet head is determined by momentum balance across 
the shock front. The shock has a Lorentz factor ~ T in the lab, and in the shock frame 
the kinetic pressure (magnetic, or radiation) p' ~ e' oc L/(r 2 r 2 ) on one side must balance 
the ram pressure from the external matter on the other side, which in the shock frame 
is PextT 2 . Their equality defines the bulk Lorentz factor of the shock, T oc r -1 / 2 , which 
is now decelerating. This continues as long as the jet continues being fed at constant 
luminosity, for t < t ac . When t = t ac is reached, independently of the initial Lorentz 
factor value, the jet will have decelerated to a value T d , reached at a deceleration radius 
r d given by 



At this time t ac ~ r d /cT d in the source frame, after feeding of the jet ceases (at this 
radius, the frozen-in pair density in the jet still exceeds the nominal Goldreich- Julian 
critical density, so the MHD assumption would remain valid). We have assumed a typical 
z ~ 20 minihalo or protogalaxy gas density of n ext ~ 100 cm -3 , e.g. Madau and Rees 
(2001). At the radius r d this density is of order of or higher than that of a possible 
& 10 _3 M Q /yr stellar wind. Beyond this deceleration radius, for a constant external 
density the jet decelerates in the energy conserving regime (e.g. Blandford and McKee, 
1976) at the self-similar rate T oc r -3 / 2 . 




cm, 



(6) 
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4. Radiation properties 

Annihilation photons .- The annihilation photons escaping from the pair photosphere at 
r a given by eq. (4) appear, in the observer frame, with a peak energy of 

Kn * T 7^\ ^ 50 keVL;/>- 8 1/ 5 2 r; /4 (20/[l + z]) (annihilation). (7) 
[l + z) 

Inverse Compton scattering in such photospheres will generally lead also to a high energy 
power law extending as N(E) oc E~ 2 above the peak (Pe'er et al. 2006). In our case, 
however, there will also be up-scattering of annihilation photons in the drag region 
r a & r ^ r p . This will depend on the scattering optical depth, which is of order unity 
just below r a and is very small above r a , but increases to t± ~ 1 at r < r p where pair 
formation sets in. One can expect a significant component of upscattered photons from 
the two radii where t± ~ 1, namely from r a at energies 0.12r a m e c 2 /(l + z), and from r p 
at energies ~ T p m e c 2 /(1 + z) in the observer frame, 

Kn,sc,a~ 50keVL^ 4 7 ri / 5 2 ry 4 (20/[l + z]) (from r a ) 

Kn,sc, P ~ 25 GeV Lj/ 4 7 r- 8 1/ 5 2 r; /4 (20/[l + z]) (from) r p ). (8) 

These components would appear as two humps at these energies, and would have a 
comparable energy, which is a significant fraction of the jet energy. Above r p , if the jet 
continues to accelerate for a while before slowing down, both the jet annihilation photons 
and external photons from the stellar envelope lead to pair formation, and some fraction 
of the jet energy could conceivably continue going into high energy photons, but the 
dynamics is dependent on the stellar model's dynamical behavior and photon input 
during the collapse and the jet propagation. A theoretical upper limit for the photon 
energy would be & Tim e c 2 , where is the inertial limit ~ 10 9 discussed below eq. (4), 
which however is unlikely to be reached, either because the field develops a longitudinal 
component, or because the flow may acquire some baryons (and associated electrons) 
whose inertia eventually becomes important. 

Interaction with external photons - Pairs and photons in the jet can also interact with 
external photons from the progenitor star or the accretion disk can exert an additional 
drag as well as resulting in additional spectral components. The latter could contribute 
significantly to the total spectrum if the jet Lorentz factor reaches values of at least 
~ T p ~ 10 6 . This is because the maximum boost in photon energy from the interaction 
is of order & T 2 , while the stellar or disk photon luminosity may be of order the photon 
Eddington value LecL,* ~ 10 41 M3 erg s _1 , giving a component of luminosity L ext & 
L Ed T 2 ^ 10 53 erg s -1 ^ L which could reach a substantial fraction of the jet Poynting 
luminosity L. 
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The process is complicated by the fact that the external (stellar or disk) radiation 
field will be inhomogeneous across the jet cross section, as well as depending on height. 
With the nominal values used here, for a jet opening angle 9j ~ lO^ 1 the lab-frame 
transverse Thomson optical depth of frozen-in jet pairs is tt,± & 1 at r ~ r p , so a 
drag (and boost) may apply over the entire jet cross section. However pair formation 
with external photons (and annihilation photons) above r p could introduce significant 
angle dependent optical depth effects. Neglecting such inhomogeneities, a discussion of 
an upscattered component and a simple example of a one-zone pair cascade spectral 
component at r p were discussed in MR97 for a Poynting jet from normal (pop. II) 
stars. Based on the expressions in that paper, the same processes would result, at 
r p in the present physical system, in a component from upscattered stellar photons at 

, l_/2ll/2 

Ef c ~ 250 TeV L^r^T 1 / \20/[l + z\), which would be absorbed by 77 interactions 
against intergalactic IR photons; and in a pair cascade component emerging at E°^ sc ~ 
2.5 keV L^^ft 8 5 2 ry 4 (20/[l + z]), whose luminosity would be given by a T p boost. 
However, transverse inhomogeneity effects as well uncertainties concerning the height 
dependence require detailed (and model dependent) calculations, making it difficult to 
say anything beyond the above semi- quantitative comments. 

Internal dissipation radiation - In this model we do not expect radiation from internal 
shocks because in a magnetically dominated outflow these do not arise (just as they do 
not play a significant role in the Crab wind); and if there aren't reversals in the field 
there shouldn't be internal dissipation (§3). 

External blast wave radiation - The contribution from this component is subject to fewer 
uncertainties than the previous ones. The forward shock from the eject plowing into the 
external medium produces a luminosity peaking at the deceleration radius r d where the 
shock Lorentz factor is ~ T d ~ 130 (eq.[6]). Following the usual treatment of external 
shocks, we estimate that for an external density in the VMS host environment of n ext ~ 
10 2 n 2 cm -3 the shocked external gas may build up turbulent magnetic fields to some 
fraction of the equipartition value with the post-shock thermal energy, resulting in a 
comoving field in the shocked gas of B' ~ (eB&KnextT^mpC 2 ) 1 / 2 ~ 8x 10 1 e^^n^T ^2.1 G. 
The comoving random Lorentz factors of the electrons in the shocked gas peak have a 
minimum (peak) Lorentz factor 7^ ~ e e Y d {m p /m e ) ~ 2.4 x 10 4 e e; _ir d 2.i, with a power 
law N(j') oc 7'~ p above that due to Fermi acceleration. The comoving synchrotron peak 
will be at an energy x' s y ~ (3/4)^7'^ ~ 8 x W^ A e 1 ^ 2 _ 1 n^ 2 T d21 el _ l , which in the lab 
frame is x sy ~ lO^e^n^r 4 2 ^ ^ ~ 50 keV, where r d>2 .i = L^t'^n' 1 ^ (eq.[6]). 
Thus, in the observer frame this is 

Ef y ~ 2.5 keV 6^ 1 e^ 1 Lj/ 2 7 t a J. 5 (20/[l + z]) (synchrotron), (9) 

independent of the assumed external (host) density. This synchrotron peak photon 
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energy is accidentally, similar to that for the external cascade at a nominal radius r p . 
As mentioned there, however, the cascade photon energy will be smeared by integration 
over r, whereas the external shock synchrotron energy (9) is generally fairly well defined, 
as observations of normal GRB indicate. The luminosity of this blast wave synchrotron 
component would be a substantial fraction of the Poynting luminosity (eq. [3]), and it 
would peak at the deceleration time t d ~ t a (eqs.[6], [1]) redshifted to the observer frame, 

^ pfc ~10 5 a:^ /2 M 3 - 1/2 ([l + ^]/20) s. (10) 

This would also be the order of the time delay between the onset of the annihilation or 
the scattering/cascade components and the blast wave synchrotron component. 

In addition, the blast wave could also have a synchrotron self-Compton (SSC) com- 
ponent, from up-scattering of the synchrotron photons by the same electrons of comoving 
Lorentz factor 7' which produced them. The synchrotron peak photons x' sy (and those 
in the power law above it) would scatter in the Klein-Nishina regime, since in the rest 
frame of the electrons x' sy ^' m > 16. Thus, the upscattered comoving photons would have 
a comoving frame peak energy ~ 7^ ~ 12 GeVe e -1^2.1, and in the observer frame 
these would appear at 

^~75GeVe e ,_ 1 I^ 2 . 1 (20/[l + z]) (SSC) (11) 

This SSC component would have a typical time lag relative to the peak synchrotron 
emission, lagging by 

tf scMg ~ {r d /cT 2 d )(l + z)~Ax 10 2 t a , 3 . 5 ([l + z]/20) s , (12) 
behind the synchrotron peak (10). 

Reverse shock radiation?.- If the jet material, beyond some radius, becomes baryon- 
loaded one would expect a reverse shock. In this case, during the period when the 
external shock Lorentz factor decreases V oc r^ 1 / 2 before reaching the value T d at the 
deceleration radius r d (eq. [6]), the reverse shock Lorentz factor would increase as Y r oc 
r 1 / 2 and could become relativistic. The reverse shock radiation would peak at r d at 
the same time t° b eak as the forward shock synchrotron component, and it should have a 
comparable bolometric luminosity. To calculate the detailed properties of this reverse 
shock would, however, require knowledge about the baryon contamination level of the 
magnetic ejecta, which is highly speculative. Nonetheless, if the properties of the rare 
prompt optical flashes of lower redshift GRBs may be extrapolated to population III 
redshifts, we can get a speculative estimate. Taking as an example either GRB 0803 19B 
at z — 0.937 with m v ~ 5 or GRB 050904 at z — 6.29 with m v ~ 10, taking the ratios 
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of the squares of the luminosity distances one would expect an optical/IR flash of up to 
~ 13 magnitudes, on similar timescales as the above X-ray and GeV components. On 
the other hand, if the baryon contamination was very low, a reverse shock would not 
form, since the Alfvenic sound speed would approach the speed of light (e.g. Giannios 
et al. (2008)). 

Afterglow.- After the blast wave radiation has peaked, in an approximately uniform 
external density the Lorentz factor would decrease as T oc r~ 3 / 2 and the usual afterglow 
would set in, with a power law time decay of the flux and a softening of the spectrum. 
The properties of this afterglow would be fairly standard, except for any differences in 
the absorption properties of the metal-poor host environment. 

5. Discussion 

In this paper we have explored in some detail the spectral and temporal properties of 
high redshift population III GRBs within the context of a Poynting dominated relativistic 
jet model. The core collapse of a pop. Ill very massive star of 250 ^ (M*/M ) ^ 10 3 
will result in an intermediate mass black hole and a temporary accretion torus, and for 
fast rotating objects the extraction of rotational energy from the black hole could power 
a Poynting dominated outflow. The largest uncertainty, in such a model, is the value of 
the Poynting luminosity L and its time-dependence. The default assumption made here 
is that, for a constant magnetic a viscosity in the torus, L is approximately constant for 
1/a times the free-fall time from the boundary of the star. This is of course uncertain 
because the efficiency of field build-up is unknown, as is whether the magnetic stresses 
get up to a significant fraction of the gas pressure P. The quantity that determines the 
jet luminosity is the field strength around the hole, which depends on the peak density 
(and peak pressure P) near the inner part of the disk. In the standard alpha model 
the density in the disk goes as r~ 3//2 . However, for a radiation dominated gas (which is 
more compressible, even if the radiation is trapped on the relevant timescale) it would in 
principle be possible for the density law to be closer to r~ 3 (and this could happen if the 
effective a were to decrease towards small r). Thus, the disk could have a higher peak 
density (and steeper profile) whatever the initial stellar profile was - and the jet could 
have a much higher luminosity for a shorter period. Our timescale estimates are also 
subject to uncertainty. We took nominally the star to be rotating at half the break-up 
speed, but the stellar angular velocity could be non-uniform. If the outer regions were 
rotating more slowly than we assumed, the the disk would obviously be smaller, but 
the times scale would still be the free-fall time, so the accretion time could be up to a 
factor a = lO" 1 ^.! shorter. Thus, there are large uncertainties in the timescales and 
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luminosities we have estimated, which could be off in either direction. Nonetheless, the 
simplifying assumption made here may be appropriate considering the preliminary state 
of knowledge about population III stars. 

The spectrum of the 'prompt' emission within the first day is shown to extend 
from soft X-rays to the multi-GeV range, with a characteristic time evolution. As a 
rough estimate, the luminosity of the annihilation photosphere (mainly in the form of 
photons) can be taken to be of the order of the remaining kinetic luminosity L/2, with 
an efficiency of conversion into annihilation photons of r] 1 ~ 1/2. The annihilation 
component of eq.(7) peaks at 50 keV and extends as N(E) oc E~ 2 up to & 3km e c 2 in 
the comoving frame, or ~ 1 MeV in the observer frame. The the spectral energy flux per 
decade E 2 (dN/dE) is the total energy times l/(\n(E max /E min ) = l/m(20) ~ 0.31. The 
radiation falling in the BAT band, 50-150 keV has 1/ ln(150/50) ~ 0.9, so we can take the 
spectral efficiency in the BAT range, say 50-150 keV, as t] B at ~ 0.5 x 0.3 x 0.9 ~ 0.13 ~ 
10 _1 . The X-ray flux might be even larger than this, adding a roughly comparable 
contribution from the external shock synchrotron component of eq. (9). With this 
efficiency, the predicted X-ray flux would be detectable in X-rays and hard X-rays by 
instruments such as the BAT detector on Swift, as estimated in §2. This would be 
detectable also in the GBM detector on Fermi, whose sensitivity is slightly better than 
Swift's. The GeV range Large Area detector (LAT) on Fermi has a fluence sensitivity 
for times t ^ 3 x 10 4 s of ~ 3 x 10~ 9 t 1 / 2 even with significantly less than 10% of the 
luminosity in the GeV band at t ~ 10 5 s, this component would be detectable by the 
LAT. To detect them, however, it may be necessary to adjust the flux trigger algorithms 
to respond to a low level, very extended increase in the flux. 

The spectral signature would have an initial hard, ~ 50 keV X-ray rise from the 
annihilation photons (equation [7]) lasting for about a day, with a possible extension out 
to ~ 25 GeV from up-scattering in the pair photosphere (equation [8]). There could be a 
cascade component from external photons leading almost simultaneously to soft X-rays 
in the few keV range, which is subject to considerable uncertainties. These would be 
followed, after a delay of hours up to a day, by an external shock synchrotron component 
in the few keV range (equation [9]). An inverse Compton component at energies in the 
70 GeV range (equation [11]) may also be expected, lagging by about ten minutes after 
the keV range external shock synchrotron component. If the jets acquire a non-negligible 
baryon load at some stage before the external shock, a reverse shock may result in an 
infrared flash of ^ 13th magnitude. An afterglow similar to that of lower redshift GRBs 
would follow over the next days, gradually shifting into the optical, infrared and radio 
frequency bands. 

The detection of such very high redshift GRBs would be of great value, as it might 
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be the first and perhaps the only way to trace the formation of the first generation of 
stellar objects in the Universe. It could give important information about the redshift 
at which the initial objects form, the rate at which they form, and the input of radiation 
into the Universe at those early epochs and its contribution to the reionization of the 
intergalactic medium. 

This research was supported by NASA NNX08AL40G and NSF PHY-0757155. We 
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